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Abstract 

We present the critical exponents 1^12, ^L2 (^nd for an m-axial Lifshitz point 
at second order in an expansion. We introduced a constraint involving the 
loop momenta along the m-dimensional subspace in order to perform two- and 
three-loop integrals. The results are valid in the range < m < d. The case 
m = corresponds to the usual Ising-like critical behavior. 
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Lifshitz multicritical points appear at the confluence of a disordered pliase, a uniformly 
ordered pliase and a modulated ordered pliase |1|, 0. The spatially modulated phase is 
characterized by a fixed equilibrium wave vector k^. In this phase, ko goes continuously 
to zero as the system approaches the Lifshitz point. If this wavevector has m-components, 
the critical system under consideration presents an m-fold Lifshitz critical behavior. This 
sort of critical behavior is present in a variety of real physical systems including high- Tc 
superconductors ferroelectric liquid crystals P, 0, magnetic compounds and alloys 

jSHlOi , among others. 

In magnetic systems ||TT|, the m-fold Lifshitz point can be described by a spin-^ Ising 



model on a ci-dimensional lattice with nearest-neighbor ferromagnetic interactions as well 
as next-nearest-neighbor competing antiferromagnetic couplings along m directions. This 
system can be described in a field-theoretic setting using a modified 0^ theory with higher 
order derivative terms, which arises as an effect of the competition along the m-directions. 
The Lifshitz universality class is defined by the parameters (A^, d, m), where is the number 
of components of the order parameter, d is the space dimension of the system, and m is the 
number of competing directions. 

Other examples of field theories containing higher derivative terms have been investigated 
in different physical scenarios. In cosmology, the recently proposed model known as "fc- 
infiation" describes inflation driven by higher order kinetic terms for the infiaton scalar 



field [0. Another instance which arises in quantum field theory in curved spacetime, is 



the quantization of scalar fields with high frequency dispersion relation around a classical 



gravitational background |jT3|. In that case, the higher order term accounts for deviations 
from Lorentz invariance. The modified dispersion relation might arise from an unspecified 
modification of the short distance structure of spacetime. A further generalization of this 
idea is to modify the large distance structure of spacetime allowing higher derivative terms. 



breaking Lorentz invariance in the infrared regime as well [|1^]. Thus a better comprehension 
of how to calculate arbitrary loop corrections for the Lifshitz critical behavior should give a 
clue about the proper perturbative treatment needed for a general higher order field theory. 



In this work we generalize the method recently developed for the m = 1 case |T5|,|T 
to calculate the critical exponents r]L2, and 'Jl using renormalization group techniques 
and the e^-expansion up to 0(e|), where ei = 4+^ — dis the expansion parameter in 



the perturbative analysis. We recover the results for the m = 1 case obtained in |]T6[ and 
show for the first time that the Lifshitz critical behavior reduces to the Ising-like one for 
m = 0. Thus, the Ising-like universality class (A^, d) is contained in a nontrivial way into 
the Lifshitz's {N,d,m). 

We start with the bare Lagrangian associated with the Lifshitz critical behavior. It can 
be written as a modified 0^ field theory expressed in the following form: 

L=^\VI 00 I' + ^1 V(d-m) 00 I' + 5o^| Vm 00 |' + ^^000 + ^^00' (l) 

The quartic dependence on the momenta along the m-directions will be manifest in the 
free propagator. Here we will consider the system at the Lifshitz critical point, defined by 
the values 60 = = 0. In order to compute the critical exponents, we need to calculate 



some Feynman diagrams, namely h, h, and [jT6|,|T^. Setting to = 5o = 0, 



[((A; + K')2)2 + (g + P)2]((A;2)2 + g2) 

is the one-loop integral contributing to the four-point function. 



(2) 



(q! + {klf) {ql + {klf) [{qi + q2+ p? + {{k, + fc^ + k'ff] 

is the two-loop "sunset" Feynman diagram of the two-point function, 

dd-^q^d'^-^q2d'^kid'^k2 



(3) 



/ 



{ql + {kl?) {{P - qiY + {{K' - k{)^f) {ql + {klf) 

X ^ 2 • (4) 



{qi-q2 + Pz? + {{ki-k2 + k'^^ 



is one of the two- loop graphs which will contribute to the fixed-point, and 

d'^^"'qid'^~"'q2d'^-"'q3d"'kid"'k2d"'k3 



{ql + {kl?) {ql + {klY) {ql + {k^f) [(gi + - P? + ((^i + A:^ - k'f Y] 

(5) 



{qi + q^ - pf + {{ki + k^ - k') 
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is the three-loop diagram contributing to the two-point vertex function. We then choose a 
special symmetry point in order to simplify the integrals. We set the external momenta at 
the quartic directions equal to zero, i.e. k' = k[ = ki^ = k'^ = 0, and K' = k'^ + ki^. In addition, 
for the four-point vertex, the external momenta along the quadratic directions are chosen as 
Pi-Pj ~ ^i^^ij ~ -'-); where Pi,P2,P3 are the independent external momenta, and P = Pi+P2- 
We fix the momentum scale of the two-point function through p"^ = = 1. We shall use 
normalization conditions for the vertex functions along with dimensional regularization for 
the calculation of the Feynman diagrams. 

Let us find out the one- loop integral l2- With our choice of the symmetry point, and 
introducing two Schwinger's parameters we obtain for I2: 



y ((Fp + (, + P)2)((P)2 + ,2) = I I da^da2^J d k exp(-(a, + a2){k 

X J d'^'^'q exp(-(ai + aa)?^ - 2a2q.P - a2P^). (6) 

The q integral can be performed to give 

J d'^-'^q exp(-(ai + aa)?^ - 2a2q.P - 02^^) 

= -Sd-nX{ „ )(Qi + Q2) 2 exp( ■ ) . (7) 

z z ai + a2 

For the k integral we perform the change of variables r"^ = kf + ... + k"^. Now take z = r^. 

The integral turns out to be: 

J d'^k exp(-(«i + a2){k^f) = (^5^)r(^)(ai + ca)"^. (8) 

Using Eqs. and (|^), I2 reads 

X / / daida2 exp{ ) (ai + a2)~^2~TJ . (9) 

Jo Jo ai + a2 

The remaining parametric integrals can be solved by a change of variables followed by a 

rescaling |T^. The integral is proportional to (P^) 2". Now we can set P^ = = 1. Using 

the identity 



T{a + bx) = r(a) [l + bxij{a) + Oix"^)], (10) 

where ip{z) = -j^lnT^z), one is able to perform the ei-expansion when the Gamma functions 
have non integer arguments. Altogether, the final result for I2 is: 
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^SmSd-m^i'^ ^ )r( ^ ) 



-(l + NmeL) , (11) 



where [i2]m = 1 + |('0(1)~~'^(2 — -j))' From now on, we shall absorb the factor inside 
the brackets in Eq. (0) in the definition of the coupling constant Then the redefined 
integral is: 

l2 = -(l + [t2]meL) . (12) 



Now we shall discuss the two- and three-loop integrals. We introduce a constraint among 
the loop momenta in different subdiagrams, along the quartic directions only [jl6[. We wish 
to highlight this approximation here by calculating the integral I4 for m 7^ 8. 

After our choice for the external momenta along the quartic directions, we can write I4 
in the following way: 

J {ql + {klf ) [{qi - 92 + V^f + ((^1 + k^W] ' ^ ' 

where we changed variables from /c2 —k2- We integrate first over the subdiagram 52, ^2- In 
order to integrate over k2 we introduce a constraint relating ki to k2 inside this subdiagram, 
i.e., ki is fixed into the second integral in Eq. (pISl). If the relation between the two loop 
momenta is of the form ki = —a k2 we can solve the integral in terms of a product of Gamma 
functions and a Hypergeometric function. The value a = 2 is singled out when we demand 
that the integral is given in terms of Gamma functions only. This is a natural generalization 
of the m = 1 case |T^. Using Schwinger's parameterization and setting ki = —2k2 in the 
second integral in Eq. ([I3|) we find 



Performing the integral over we obtain 

d'^-'^qi 



(14) 



dz 



iql-2zP.q, + zP^f-- [(gi+ps)'] 
Using a Feynman parameter the integral turns out to be 



(15) 
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X dz 
Jo 



yz{l - yz)P'^ + y{l - y)pl - 2yz{\ - y)p3.P 



(16) 



The integral over y is singular at y = 1 when = 0. We only need to replace the value 
y = 1 inside the integral over z [1^,17], and integrate over y afterwards, obtaining 



^4 = ^(l + 3 [i2\meL)- 



2eL 

The integrals R and li can be solved using a similar reasoning. They are given by 



(17) 



I'. 



1 1 



8 — m ei 



1 + [l2]m + 



^ 2 



(18) 



3(2 - f ) ei 



1 + 2 [i2]m + 



2 - 



(19) 



Note that the leading singularities for /2, /4 are the same as their analogous integrals 
in the pure 0^ theory. However, Jg and Jg do not have the same leading singularities for 
they include a factor of 1^2-^) ■ then introduce a weight factor for /g and /g, namely 
(1 — ^), so that they have the same leading singularities as in the pure 0^ theory. This 
has the advantage of allowing a smooth transition to the Ising-like case (m = 0) from the 
general Lifshitz anisotropic critical behavior (m 7^ 8) as we shall see next. 

The fixed point at two-loop level is given by: 



u 



6 

8 + N 



eii 1 + 



/ 4(5iV + 22) 
I [8 + NY 



- 1 M 



{2 + N) 
(8 + A^)2 



(20) 



With this fixed-point one readily obtains the critical exponents r]L2 and iyL2- 
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1 2 2 + N 
2^^ (8 + Nf 



, (2 + iV) 

+ , T\T\2 



{8 + Ny 



/ 4(5A^ + 22) 
1^ (8 + iV)2 



^2 m + 



— m 



2 + N 

8 + Ny 



(21) 



_ 1 1 2 + 
^^^-2 + 4'^8T^ 



+ 



1 (2 + A^) 
8 (8 + iV)3 



2(14iV + 40) [^2]«^ - 2(2 + N) + {8 + N){3 + N) 



el- 



Using the scaling law 7^ = 2^^2(2 — VL2), the exponent is 



, 1 2 + N 
^^ = ^ + r^8TiV 



(22) 



12 + 8A^ + AT^ + 4 [i2U (20 + TAT) 



(23) 



4 (8 + N)^ 

It should be emphasized that [i2]m is a universal amount, for the dependence on m is 
encoded in such quantity. The parameter m only appears in a explicit way at the ^(e^,) 
contribution to the index 77^2- To our knowledge, the explicit dependence on m is obtained for 
the first time at 0(ei) for r]L2- When setting (m = 1) in the formulae above, we recover the 
exponents previously reported in reference |TB|. As discussed there, the two-loop calculation 
(A^ = 1) in three dimensions yields 7^ = 1.45, in a nice agreement with the numerical Monte 
Carlo simulation = 1.4±0.06. 

The amazing fact obtained using the method outlined here is that the critical exponents 
reduce to the Ising-like ones when m = 0, for eL^e = 4 — d. This means that the 
universality class for the m— fold Lifshitz point includes the Ising-like one for this particular 
value of m in a nontrivial way. This provides a unified description of the anisotropic Lifshitz 
critical behavior (m 7^ 8, d ^ m). This is the first time that an isotropic behavior (m = 0) 
can be recovered from the most general anisotropic Lifshitz criticality. 



Note that our result for the exponent r]L2 is in agreement with Mukamel's [|T^] at 0{eD 
and is independent on m at this order. It should not be surprising that the approach fails to 
describe the d = m = 8 case, for the exponent 77^2 is divergent as can be seen from Eg. (pH) . 
The approximation made is not suitable for general isotropic cases d = m 7^ 8 as well, since 
there is no preferred directions any longer. Another treatment should be employed to obtain 
information along the m-dimensional competition axes, since the symmetry point used here 
is not suitable to find out quantities along the competing directions. 

All the results in this work follow from expanding the theory around its upper critical 
dimension. The constraint introduced along the m-dimensional subspace is equivalent to 
expand around the theory without competition, with m kept fixed. A different field-theoretic 
method has been proposed, based on the expansion around the number of the m competing 
directions pO|-p^. The m = 2,6 reported cases are in disagreement with our results. This 



suggests that the two approaches are inequivalent. 

To conclude, we have calculated the critical exponents associated to correlations along 
the {d — m)-directions perpendicular to the competition axes. This was possible because we 
introduced a constraint between the quartic loop momenta appearing in different subdia- 
grams in higher-loop Feynman graphs. The Lifshitz universality class turns out to reduce to 
the Ising-like one for the value m = at least up to the loop order considered in this work. 
In principle, the technique can be readily generalized to analyze general anisotropic Lifshitz 
type critical behavior with arbitrary powers of the Laplacian at the competing directions. 
The study of the tricritical Lifshitz points using this formalism is also worthwhile. 
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